New limit on the lepton-flavour violating decay /i"*" — )■ e"'"7 
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We present a new result based on an analysis of the data collected by the MEG detector at the 
Paul Scherrer Institut in 2009 and 2010, in search of the lepton flavour violating decay ^'^ — !> 6^7. 
The likelihood analysis of the combined data sample, which corresponds to a total of 1.8 x lO'^* 
muon decays, gives a 90% C.L. upper limit of 2.4 x lO"^'^ on the branching ratio of the ^"'" -^ e'^-y 
decay, constituting the most stringent limit on the existence of this decay to date. 
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The lepton flavour violating (LFV) decay /x — ;> 67 is 
forbidden within the standard model of elementary parti- 
cles (SM). Even with the introduction of neutrino masses 
and mixing SM predicts an immeasurably small branch- 
ing ratio {B < 10^^^) for this decay. Conversely new 
physics scenarios beyond SM, such as supersymmetric 
grand unified theories or theories with extra dimensions, 
predict branching ratios in the 10~^^ to lO"^** range il^- 
|3|. This is close to the present limit set by the MEGA 
experiment [^, B < 1.2 x 10~^^, which places one of 
the most stringent constraints on the formulation of such 
theories. Observation of /i — > 67 therefore would be an 
unambiguous signature of new physics, while improve- 
ments on the existing limit would stringently constrain 
many of the new physics scenarios beyond SM. 

The MEG experiment [5|, l6( covers a 10% solid angle, 
centred around a thin muon stopping target (205 /im- 
thick polyethylene) and is composed of a positron spec- 
trometer and a photon detector in search of back-to-back, 
monoenergetic, time coincident photons and positrons 
from the two-body /i^ — > e"'"7 decays. The positron spec- 
trometer consists of a set of drift chambers (DC) 7] and 



scintillation timing counters (TC) [8j located inside a su- 
perconducting solenoid with a gradient field [9j along the 
beam axis, ranging from 1.27 Tesla at the centre to 0.49 
Tesla at either end. The photon detector [10|, located 
outside of the solenoid, is a homogeneous volume (900 £) 
of liquid xenon (LXe) viewed by 846 UV-sensitive photo- 
multiplier tubes (PMTs) submerged in the liquid. The 
spectrometer measures the positron momentum vector 
and timing, while the LXe detector is used to reconstruct 
the 7— ray energy as well as the position and time of its 
first interaction in LXe. All the signals are individually 
digitized by in-house designed waveform digitizers based 
on the multi-GHz domino ring sampler chip (DRS) [llj . 
The PSI 7rE5 beam line is used to stop 3 x 10^ posi- 
tive muons per second in the target. The residual polar- 
ization of the decaying muons along the beam axis was 
measured to be (P) = —0.89 ± 0.04. The background to 
/i+ — > e"'"7 decay comes either from radiative muon de- 
cays /i+ — >■ e+i^P7 (RMD) in which the neutrinos carry 
away little energy or from an accidental coincidence of 
an energetic positron from a normal Michel decay with a 
7— ray coming from RMD, Bremsstrahlung or positron 



annihilation-in-flight. The accidental coincidences are 
the dominant background in this experiment. 

The MEG detector response, resolutions and stabil- 
ity are constantly monitored and calibrated. The pho- 
tomultipliers (PMTs) of the LXe detector are calibrated 
daily by LEDs and a-sources immersed in the liquid 12 1 . 
The energy scale and resolutions of the LXe detector are 
measured over the energy range of 4.43 to 129.4 MeV 
using 7— rays from a radioactive Am/Be source, (p, 7)- 
reaction u sing a dedicated Cockcroft- Walton accelera- 
tor (CW) [l3|, and n'p charge exchange and radiative 
capture reactions (CEX). A 9MeV-7 line from the cap- 
ture in nickel of neutrons from a pulsed and triggerable 
deuteron-deuteron neutron generator allows one to check 
the stability of the LXe detector even during data-taking. 
The relative time between the TC and LXe detector is 
monitored using RMD and 27-events from ^lB{p, 27)^|C 
reactions. 

The fi~^ — >■ e^7 trigger requires the presence of a high 
energy 7— ray in the LXe detector and a hit on the timing 
counters within a 20 ns window together with an approx- 
imate back-to-back topology. Pre-scaled monitoring and 
calibration triggers are also recorded. A more detailed 
description of the MEG detector can be found in Ref. [6J . 

The results presented in this paper are based on data 
collected in 2009 and 2010 (for a total of 1.8 x lO" //+- 
decays in the target); the 2010 statistics are about twice 
that of 2009. All sub-detectors were running stably dur- 
ing these periods. The 2008 data [Q] are not used in 
this analysis because of their limited statistics and de- 
tector performance. In 2010 a DRS upgrade resulted in 
an improvement in the time resolution while an increase 
in noise in the DC, due to a deterioration of the HV 
power supplies, caused slightly worse positron tracking 
resolutions. 

We adopted a likelihood analysis method combined 
with a blind procedure on examining the data: events 
close to the signal region were kept hidden (blind region) 
until all the analysis procedures had been completely de- 
fined. The probability density functions (PDFs) needed 
for the likelihood analysis were constructed using the 
events outside of the blind region (side-bands) . 

Several improvements to the analysis have been intro- 
duced since the presentation of the preliminary result 
based on the 2009 data lj| and also implemented in the 
current 2009, as well as in our 2010 analyses. These im- 
provements include a new alignment technique for the 
DC system; an improved experimental evaluation of the 
spectrometer performances; a better understanding of 
the gradient magnetic field; improvement in the relative 
alignment of the photon detector and the positron spec- 
trometer by means of cosmic ray muons; adoption of a 
more commonly used statistical method (profile likeli- 
hood); a constraint on the background rates in the like- 
lihood analysis from the data in the side-bands. 

The kinematic variables used to identify the fj,~^ — >■ e^7 



decays are the 7-ray and e+ energies {E^, Ee), their rel- 
ative directions {Oe-y, (pe-y) UM ^^'^ emission time (^67)- 
The offline event selection requires at least one e^-track 
reconstructed in the spectrometer and pointing to the 
target, with minimal quality cuts applied. The blind re- 
gion is defined by A8 < E^ < 58 MeV and |ie^| < 1 ns. 

The positron track reconstruction in the spectrometer 
is based on a Kalman filter technique 16|. Effects of 
multiple scattering and energy loss in the detector mate- 
rials in the presence of the non- uniform magnetic field are 
taken into account. Internal alignment of the DC is ob- 
tained by tracking cosmic ray muons without a magnetic 
field and by minimizing the measured residuals in a man- 
ner independent of the initially assumed alignment 17 1 . 
The absolute position of the DC system is based on an 
optical survey. 

The magnetic field of the spectrometer was measured 
at the beginning of the experiment and only its major 
component along the beam axis is used in the analysis to 
avoid possible misalignment errors from the Hall probes; 
the other minor components are deduced from the deriva- 
tives of the measured primary component using Maxwell 
equations together with boundary conditions at a sym- 
metry plane at the magnet centre where the minor com- 
ponents are nearly zero. This magnetic field map agrees 
to within 0.2% with the field computed for the geometry 
and currents of the spectrometer coils. 

The resolutions of the positron track direction are es- 
timated by exploiting tracks with two full turns in the 
DC. Each turn is treated as an independent track and 
the resolutions are extracted from the difference between 
the two reconstructed sections. The energy resolution is 
evaluated by fitting the kinematic edge of the Michel de- 
cays and is well described by a sum of three Gaussians 
with resolutions of 0.31 MeV, 1.1 MeV and 2.0 MeV for 
the core (80%) and the two tail (13% and 7%) compo- 
nents, respectively for 2009 and 0.32 MeV, 1.0 MeV and 
2.0 MeV for the core (79%) and the two tail (14% and 
7%) components, respectively for 2010. 

The decay vertex coordinates and the positron direc- 
tion at the vertex are determined by extrapolating the 
reconstructed track back to the target. The 7-ray direc- 
tion is defined by the line connecting the decay vertex to 
the 7-ray conversion point measured by the LXe detector. 

A geometrical correlation exists between errors on 4>e 
at the vertex position and E,, , which is measured by using 
the two-turn method and is perfectly reproduced by the 
Monte Carlo (MC) simulation. The </>£ -resolution has a 
^e-dependence due to the correlation and has a minimum 
at (f>e = 0, where it is measured to be a^^ = 6.7(7.2) mrad 
for 2009 (2010) data ^. The 6'c-resolution is measured by 



^ From here on we will quote in parentheses the value in the 2010 
data when different from that in 2009. 



the two-turn method to be ug^ = 9.4(11.0) mrad. The 
resolution on the decay vertex coordinates is also deter- 
mined by the two-turn method; along the beam axis it 
is described by a Gaussian with ct^ = 1.5(2.0) mm while 
in the vertical direction it is described by the sum of two 
Gaussians with ay = 1.1mm for the core (87%(85%)) 
and ay = 5.3(4.8) mm for the tail. 

The determination of the photon energy E-^ in the LXe 
detector is based on the sum of the number of scintilla- 
tion photons detected by the PMTs; correction factors 
take into account the different PMT geometrical accep- 
tances. Due to its geometry the detector response is 
not totally uniform over the photon entrance window; 
this is corrected for by using 7-lines from CW and CEX 
reactions. The absolute energy scale and resolution at 
the signal energy E-y = 52.8 MeV are determined by the 
GEX measurement; the resolution ur, extracted from a 
Gaussian fit to the high energy side of the spectrum, 
depends also on the depth {w) of the 7— ray conversion 
point from the photon entrance surface of the LXe de- 
tector: CTR = 1.9%(u; > 2cm) and 2.4%(m; < 2cm). The 
3D-map of the measured resolutions is incorporated into 
the PDFs for the likelihood analysis. 

The photon energy scale and the resolutions are cross- 
checked by fitting the background spectra measured in 
the side-bands with the theoretical RMD spectrum folded 
with the detector resolutions; the resolutions during the 
run are well represented by the CEX evaluations and the 
systematic uncertainty of the iJ-^-scale is estimated to be 
~ 0.3%. Since MEG operates at a high beam intensity, 
it is important to recognize and unfold pile-up photons. 
For each event the spatial and temporal distributions of 
the PMT charge are studied to identify photon pile-up 
in the LXe detector; in case of positive identification, 
corrections to the PMT charges are applied. Cosmic ray 
events are rejected using their characteristic PMT charge 
distribution. 

The position of the first interaction of the 7-ray in 
the LXe detector is derived from the light distribution 
measured by the PMTs close to the region of the energy 
deposition by fitting the distribution with the expecta- 
tion. The position resolution in the plane of the photon 
entrance window is measured to be 5 mm in a dedicated 
CEX run with a lead slit-collimator placed in front of the 
LXe detector, while the resolution along the depth w of 
6 mm and the position dependence of the resolutions are 
evaluated by a Monte Carlo simulation. 



The resolutions on the relative directions 



r) 



are derived by combining the relevant resolutions of 
positrons and photons discussed above; the results are 
14.5 (17.1) mrad for 6'e^ and 13.1 (14.0) mrad for 0e^. 
The relative time te-y is derived from the two time mea- 
surements by the LXe detector and the TC, after cor- 
recting for the length of the particle flight-path. The 
associated resolutions at the signal energy 146(122) ps 
are evaluated from the RMD peak observed in the E^ 



side-band; a small correction takes into account the E^- 
dependence measured in the CEX calibration runs. The 
position of the RMD-peak corresponding to te-y = was 
monitored constantly during the physics data-taking pe- 
riod and found to be stable to within 15 ps. 

A likelihood analysis is carried out for events in a 
portion of the blind region (analysis region) defined by 
48 < £:^ < 58 MeV, 50 < Ee < 56 MeV, |ie^| < 0.7 ns, 
\Oe-y\ < 50 mrad and |(/)e^| < 50 mrad. These intervals in 
the analysis variables are between five and twenty sig- 
mas wide to fully contain the signal events and also re- 
tain some background events. The best estimates of the 
numbers of signal, RMD and accidental background (BG) 
events in the analysis region are obtained by maximizing 
the following likelihood function: 



„-N (Wrmd-(Wrmd))^ (Wbg-("bg>)^ 
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where Xi — {-E-y, i?e,ie7j^e7 7 '^'67} is the vector of ob- 
servables for the «-th event, A^^^ig, A^rmd and TVbg are 
the expected numbers of signal, RMD and BG events, 
while S, R and B are their corresponding PDFs. N = 
iVsig + iVRMD + iVsG and 7Vobs(= 311(645)) is the ob- 
served total number of events in the analysis window. 
(A^rmd)(= 27.2(52.2)) and (iVBG)(= 270.9(610.8)) are 
the numbers of RMD and BG events extrapolated from 
the side-bands together with their uncertainties crRMD(= 
2.8(6.0)) and crBG(= 8.3(12.6)), respectively. 

The signal PDF S{xi) is the product of the PDFs for 



Ee, 9, 



e-yi 067 and tej, which are correlated variables, as 



explained above, and the E^ PDF. The PDFs properly 
incorporate the measured resolutions and correlations 



among E^ , 9e 



and te-y on an event-by-event basis. 



The RMD PDF R{x,) is the product of the same te7-PDF 
as that of the signal and the PDF of the other four cor- 
related observables, which is formed by folding the the- 
oretical spectrum with the detector response functions. 
The BG PDF B{xi) is the product of the five PDFs, each 
of which is defined by the single background spectrum, 
precisely measured in the side-bands. The dependence 
of the resolutions on the position of the 7-ray interaction 
point and on the positron tracking quality is taken into 
account in the PDFs. 

A frequentist approach with a profile likelihood-ratio 
ordering 18J, |l9| is used to compute the confidence inter- 



vals on A^sig: 



Ap(A^sig) 



CjNsig, A^RMD(A^sig), A^BG(iVsig)) 
^(A^sig,^RMD,^BG) 



where the hat and double hat denote the best estimates 
maximizing the hkehhood for floating and fixed -/Vgig, re- 
spectively. Other, independent analysis schemes based 
on averaged PDFs without event- by- event information 
or Bayesian approach were also used and found to be 
compatible with the analysis presented here to within 10 
to 20% difference in the obtained branching ratio upper 
limits. 

In order to convert TVgig into a branching ratio value 
the normalization relative to the Michel decay is com- 
puted [GJ by counting the number of Michel positrons 
passing the same analysis cuts. This is accomplished 
by means of a pre-scaled Michel positron trigger enabled 
during the physics data-taking. A correction to the pre- 
scaling factor due to positron pile-up in the TC is taken 
into account. Another method for computing the nor- 
malization uses RMD events in the Ej side-band and the 
theoretical branching ratio of the RMD. The normaliza- 
tions calculated by these two independent methods are 
in good agreement and are combined to give the normal- 
ization factor with a 7% uncertainty. 

The sensitivity of the experiment with a null signal hy- 
pothesis is evaluated by taking the median of the distri- 
bution of the upper limit on the branching ratio obtained 
over an ensemble of toy MC experiments. The rates of 
RMD and BG events, as measured in the side-bands, are 
assumed in the simulated experiments. The branching 
ratio sensitivity at 90% confidence level (C.L.) is found 
to be 3.3 X 10-12 (2.2 x lO-^^) for the 2009 (2010) data 
sample and 1.6 x 10^ ^^ when 2009 and 2010 are com- 
bined. These sensitivities are consistent with the upper 
limits obtained by the likelihood analyses in several com- 
parable analysis regions of the 1,,^ side-bands. 

After calibrations, optimization of the analysis algo- 
rithms and background studies in the side-bands are com- 
pleted, the likelihood analysis in the analysis region is 
performed. In Figures [T] we present the distributions, 
for the 2009 and 2010 data samples respectively, showing 
the events seen in the analysis region projected in the E-y 
vs Ee and te-y vs cosOe^ planes, Ge-y being the opening 
angle between the 7-ray and the positron. In plots (a) 
and (c) selections in tey and cosGg-y, each of which is 
90% efficient on the signal, are applied {\tej\ < 0.28ns 
and cos Qe-y < —0.9996) ; in plots (b) and (d) a selection 
in Ee which is 90% efficient on the signal and a selec- 
tion in Ej which is 73% efficient on the signal inside the 
analysis window are applied (52.3 < Eg < 55MeV and 
51 < ^^ < 55MeV). The contours of the signal PDF 
are also drawn and a few events with the highest signal 
likelihood are numbered in a decreasing order of relative 
signal likelihood, S/{fKR + fBB), /r = 0.1 and /b = 0.9 
being the fractions of the RMD and the BG measured in 
the side-bands, respectively. High signal likelihood events 
were thoroughly checked and found to be randomly dis- 
tributed in time and detector acceptance. 

The observed profile likelihood ratios as a function of 
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FIG. 1: Event distributions in the analysis region of (a) E.y 
vs Ee and (b) t^-y vs cos Qey for 2009 data and of (c) E^ vs 
Ee and (d) te-, vs cos Qey for 2010 data. The contours of the 
PDFs (1-, 1.64- and 2-ct) are shown, and a few events with 
the highest signal likelihood are numbered for each year. (The 
two highest signal likelihood events in 2010 data appear only 
in (c) or (d).) 



the branching ratio for 2009, 2010 and the combined data 
sample are shown in Fig. [2] [20[ . The analysis of the full 
data sample gives a 90% C.L. upper limit of 2.4 x IQ--'^^, 
which constitutes the most stringent limit on the exis- 
tence of the fj,'^ -^ e~^"f decay, superseding the previous 
limit by a factor of 5. The 90% C.L. intervals as well as 
the best estimate of the branching ratio for 2009 and 2010 
data separately are also given in Table HI The 2009 data 
set, which gives a positive best estimate for the branch- 
ing ratio, is consistent with the hypothesis B — with 
an 8% probability. 

The systematic uncertainties for the parameters of the 
PDFs and the normalization factor are taken into ac- 
count in the calculation of the confidence intervals by 
fluctuating the PDFs according to the uncertainties. The 
largest contributions to the systematic uncertainty, which 
amount to a shift of about 2% in total in the branching 
ratio upper limit, come from the uncertainties of the off- 
sets of the relative angles, the correlations in the positron 
observables and the normalization. 

The MEG experiment continues data-taking and is ex- 
pected to explore the /i+ — > e+7 decay down to a branch- 
ing ratio sensitivity of a few times lO-i'^ in the next few 
years. 
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FIG. 2: Profile likelihood ratios as a function of tfie /i^ — > e+7 
brandling ratio for 2009, 2010 and the combined 2009 + 2010 
data sample. 



TABLE I: Best fit (Bfu), lower (LL) and upper limits (UL) 
at the 90% C.L. of the branching ratio for the 2009, 2010 and 
combined 2009 + 2010 data sets. 



Data set 



Bfit 



LL 



UL 



2009 


3.2 X 10"^^ 


1.7 X 10"" 9.6 X W'^^ 


2010 


-9.9 X 10^" 


1.7x10^^2 


2009 + 2010 


-1.5 X 10"" 


2.4 X 10-^2 
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